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Abstract. A complete classification of linear differential operators possessing 
finite-dimensional invariant subspace with a basis of monomials is presented. 



1. Introduction 

One of the old-standing problems in the theory of special functions is the classi- 
fication of all linear differential operators that admit an infinite sequence of orthog- 
onal eigenvectors in the form of polynomials. See [2] for an overview. In 1929 S. 
Bochner [1] had solved this problem for second order differential operators on the 
complex (real) line. Therefore we name this problem the "Bochner problem" . 

The main purpose of the present paper is to study a more general problem, 
which we name (following [4]) the "generalized Bochner problem". We ask for a 
classification of all linear differential operators, which possess a finite-dimensional 
invariant subspace of polynomials. It turns out this problem is rather sophisticated. 
However, things get much more tangible when, instead of invariant subspace in 
polynomials, we require that the (finite-dimensional) invariant subspace has a basis 
of monomials. This problem is solved completely (Section 3). In [4] a particular 
case of this problem was solved, namely if the invariant subspace is the linear space 
of polynomials of degree not higher than some fixed integer. The classification of 
the linear operators possessing such an invariant subspace was given through the 
universal enveloping algebra of the algebra sl2 taken in a special representation by 
first-order differential operators. 

The results of Section 3 have an impact to the general problem. This is pre- 
sented in Section 4. In Section 5, we give the explicit expressions for the second- 
order differential operators T2, possessing a finite-dimensional invariant subspace 
in monomials. These operators are of great interest for finding explicit solutions to 
the Schrodinger equation 

(1.1) (^.^ + v{x))^{x)^\^{x) 

ax 

since the eigenvalue problem for the operator T2,T2(p = X(p can be reduced to the 
Schrodinger equation by a change of the variable, x' — x' {x), and introducing a new 
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(gauge-transformed) functfon, ^' = !/?exp(— o(a;)). Hence the operators described 
in Section 5 lead to a special class of quasi-exactly-solvable Schrodinger equations 
complementary those described already [4] . 

Finally we make some concluding remarks in Section 6 concerning the case where 
the powers of x are not natural numbers. 

Though in this paper we work over the complex numbers, the main results hold 
for real numbers as well. 

2. Generalities 

Here we state some general results concerning differential operators which leave a 
finitc-dimcnsional subspacc of C[a;] invariant. The results of this section can easily 
be extended to more variables. Let D denote the algebra of linear (finite-order) 
differential operators on C[a;] with polynomial coefficients. We denote the symbol 

Let y be a finite-dimensional subspace in <C[x\ of dimension n, and let T : 
C[x] C[x] be a linear operator. It is easy to show that any linear operator T can 
be written as an infinite order linear differential operator: 

oo 

This can be proved easily by an inductive construction: 
Po = T(l), P,=T{x)-xPo,... . 

It follows immediately that the action of T on F can be represented by a finite-order 

differential operator, since V C JPk (where JPfe denotes the space of polynomials of 
degree not higher than k, and k sufficiently large), and hence for Tj^: 

k 

we have that Tfe(v) = T{v) for all v eV. 
So the following proposition is immediate: 

Proposition 2.1. Let V C C[x] be a subspace of dimension n, and let Dv denote 

the algebra of differential operators that leave V invariant (i.e. Tk G =^ Tk{V) C 
V). Then Dv is isomorphic to the semi-direct product o/ End(y) and I, where 
End(y) is the algebra of linear operators on V, and the ideal I is the algebra of 
differential operators which annihilate V. □ 

We can paraphrase this proposition in the following way. As representation on 
V, Dv is the full n x n-matrix algebra, and the kernel of this representation is an 
infinite dimensional ideal inDv- 

Another consequence of this observation is. that for any element v € C[x] and 
A e C, we can find a differential operator T such that Tv = Xv. 

In this paper, we mainly discuss the case that V is graded. Now C[a:] is a 
graded algebra by putting dcg(.T'') = k. Let us take V a graded subspace. This 
means exactly that V has a basis of monomials. Hence we assume that V = 
{x^^,x^^, . . . ,a;'"), which we abbreviate toV= (x^), I = {ii,i2, ■ ■ ■ ,in}- 
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The fact that C[x] is graded has as a consequence that End(C[a;]) is also graded, 
putting deg(T) = m, if T{x^) G (a;'"*"™) for all i. Moreover for graded V it follows 
that 5)y is also graded. So, in order to describe the structure of Sy, it is sufficient 
to describe the homogeneous components of Sy. This is the concern of the next 
section. 

3. Algebras leaving a space of monomials invariant 

As before, we arc interested in finite-order linear differential operators T, such 
that T{V) C V. We assume T to be graded of degree m and order k, which means 
that 

k 
i=0 

Here c/j ^ and m G Z. Moreover Cj = for i + m < 0. In particular, we see that 
if the degree of T is negative, say dcg(r) = — m, m > 0, then the order of T is at 
least m. 

The following lemma plays a crucial role in our classification. 

Lemma 3.1. Let T be a differential operator of degree m and order k. 

(1) Suppose m > 0. Then there exist numbers ai, . . . , afc G C such that 

T = Q,.T™(x9 — ai){xd — a2) ■ ■ ■ [xd — ak). 

(2) Suppose m < 0, hence k > —m. Then there exist numbers ai, . . . , ak+m S 
C such that 

T = cud^"'-(xd - ai)(xd - ■ ■ ■ {xd - afe+m)- 
Proof. If m > it is clear that any differential operator of degree m is of the form 

k ' k 

T = x"^ ^ CiX^d^. Similarly, if m < 0, T can be put in the form T = 9^™ Yl CiX^d'^ 

j=0 2=0 

with Cfe — Cfe. So it suffices to prove the lemma for m = 0. But then we have 
T{x"') = Pk{a)x°', where Pk is a k-th order polynomial. So Pk{oi) = Ck{c( — ai){a — 
0:2) • • • (a — ak) for some at G C. Now 

{{xd - Q!i) • • • {xd - ak)){x") = {a- ai){a - 02) • • • (a - Uk), 

From this it follows that T = Ck{xd — ai){xd — a^) ■ • ■ {xd — ak), since the repre- 
sentation of !D on C[.t] is faithful. □ 

It is easy to see that this factorization is unique up to the ordering of the factors 
xd — ai. However, a certain ordering of the factors has no meaning, since these 

factors commute. 

Above representation of T given by lemma 3.1 is very convenient, when we 
study the operators T which leave V = (a;'^, . . . ,x^") invariant, i.e. T G Dv- We 
introduce the following notation. For I = {ii,i2, . ■ . , in} and m G Z, we put 

/C™) = {iel\i + m>0 and i + m ^ 7}. 

If deg(T) = m and T G then it is clear that with necessity T{x'^) = for 
i G J(™\ since T{x') = c ■ but x'+'^ ^ V, so c = 0. This leads to the 

following 
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Theorem 3.2. Let V = and T a finite-order differential operator 

such that deg(T) = m. Suppose I^"^^ = {ai, . . . , afc}. Then T e Dv, if o-nd only if 

T = f-{xd- ai){xd - 02) ■ ■ ■ {xd - ak), 

where T is some differential operator of degree m. 

Proof. The if-part is trivial. So assume that T G Dy, and suppose it has order s. 
According to lemma 3.1, for m > 0, T can be represented in the form 

T = c- x'^ixd - 0i){xd - 02) ■ ■ ■ {xd - f3s) (c e C, c ^ 0). 

We need that T{x^) — for i £ /(™). On the other hand, wc have 

T{x') = c{i - - /?2) • • • - /3.)x*+". 

Hence it follows that {ai, . . . , au} C . . . , /3s}. After rearranging the /3's we 
find 

T = c- x'^ixd - f3i)(xd - /32) • • • {xd - (3s-k){xd - ai) ■ ■ ■ {xd - ak). 

So for m > the proposition is proved; T = c- x™' {xd — Pi){xd — P2) • ■ • {xd — Ps-k)- 
For m < the proof is similar. □ 

REMARK. Prom the previous proposition it follows that the order of T e S)y 

with dcg(T) = m is at least k, where k is the number of elements in In fact, 

up to a scalar coeflScient, the element of order k is unique: 

T = x"'{xd - ai) ■ ■ ■ {xd - ak), 

where {ai, . . . , ak} = I^"^^. 

If m = 0, = 0, and we have T = 1. So all differential operators of degree 
are in these elements form a commutative subalgebra of S)y generated by xd. 

4. The case of polynomial subspaces 

In section 3, we performed the classification of differential operators with an 

invariant subspacc V that has a basis of monomials. If V has no basis of monomials, 
Dv is not be graded, but only filtered. This causes a major difficulty. However, 
considering the corresponding grading, we still can deduce some properties of 2)y 

in this case. 

So, let V have a basis of the form 

X^^ + Ciix'^-'^ + . . . , X''' + C21X'-'-'^ + . . . , ... , X'^ + Cnix'^''^ + . . . 

We can assume that all ij are different. The graded space V{g) associated to V is 
{xi\... ,x^-). 

Let T G S) of order k be of the form 

m 
i=—k 

with T^') of degree i, and T^™) ^ 0. (An operator of order k has degree —k or 
higher, the term of degree —k being a multiple of 9'^). We call T^™) the associated 
graded operator. Now it is easy to prove 
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Theorem 4.1. Let V, V{g), T andT^""^ be as above. 
IfTG Dv, then G ^v{g)- 

Proof. Suppose T e Dv, and consider T{x^^ + Cjix^^~^ + ...). If m = there is 
nothing to prove, so consider m 7^ 0. Let / = {ii,i2, ■ • • ,in}- If ij G -f^™^ then 
should be 0, since no term of T can cancel this term. Hence we find 
exactly that T(™)(a;'J) = for all ij e i^™), i.e. T^'") e □ 

Prom this we derive an easy corollary: 
IfTk G D of order k possesses an infinite number of (linearly independent) eigen- 
vectors, then T = E°=-fc ^^'^ ^ef'^ee o/T^^) = i and r(°) ^ 0. 

A similar reasoning can be performed for the part of T with minimal degree, but 
this seems to give not much information. 

5. Classification of second order differential operators 

5.1. Generic situation. Now we proceed to 2-nd order differential operators T2, 
which admit a finite-dimensional invariant space of polynomials with a basis of 
monomials. We are interested in this problem in connection with finding explicit 
solutions to the Schrodinger equation (1-1). As mentioned in the Introduction, this 
involves some transformations see [3] . These transformations are of 2 types: the first 
is a change of basis, and the second is called "gauge" transformation, which amounts 
to changing T to gTg~^, where g is a non-zero function. To delete some ambiguity 
in our spaces of monomials, we impose 2 conditions on y = {x'^^ , a;'^ , . . . , a;'" }: 

(1) We assume that 1 G F (it removes an ambiguity resulting from gauge 

transformations) . 

(2) We assume that gcd(ii,Z2, . . . ,z„) = 1, i.e. that the powers have no com- 
mon factor (it removes an ambiguity resulting from changes of variable). 

For the classification of differential operators in S)y these assumptions do not make 
much difference. liT2 = T^^^ , where the degree of T^*) is i, these two assumptions 
effect only the terms T^-^) and T^-^). 

So let us start the classification. Suppose T2 G Dy, T2 = Em^^"^ where T^™) 
has degree m and order less or equal to 2. If T^™' is non-zero, then according to 
proposition 3.2, /(™) contains 0,1 or 2 elements. We distinguish these 4 cases: 

(1) For all m > 0, 7^"*^ contains more than 2 elements. So T2 contains no 
terms of positive degree. Hence T2 preserves (l,x,a;^,... ,a;") for all n. 
This type of operators is already studied in [4] called there exactly-solvable 
operators), and we do not repeat it here. 

(2) 7^"*^ is empty. This can only happen, if to = 0. But this case is trivial, 
since all operators of degree are in Sy. Hence degree contributes to T2 
the operator aix'^d^ + a2xd + 03. This part we call trivial, and is always 
present in T2. 

(3) I^"^^ contains one element. Let us assume here that «i > ^2 >•••> in- 
Then we have ij = ij-i + m so that ij = {n — j)m (since we assumed that 
I E V , so in = 0) . But we also assumed that the ij have no common factor, 
so it follows that to = 1, and hence V = {1 Hence we are 
in the case that is extensively discussed in [4]. 
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(4) contains 2 elements, m > 0, and no /^'^ for I > contains one element. 
Suppose ii and i2 are the 2 elements in Then the set {ii,i2, ■ ■ ■ ,in} 
is of the following specific form: 

ia = ii - m, = is - m, ... , i2r-i = hr-a - m 

and 

ii = 12 — m, ie = i4 - m, . . . , i2s = «2.s-2 — m 

We call H, is, . . . , t2r-i and *2, «4, • ■ ■ , «2s chains with step m and length r 
and s, respectively. In general, with no special relation, by which we mean 
that / can be split into two chains in exactly one way, the most general 
second-order operator inDv is 

T2 = aix"^{xd — ii){xd — 12) + a2X^d^ + a^xd + 

and T2 contains an extra term in two cases: 

(a) If m = 1, (and therefore i2s = 0), we get an extra term 
a5d{xd - i2r-i)- 

(b) If m=2 and {^2^, i2r-i} = {0, 1}, we get an extra term a^d"^. 

All this can easily be proved, by examining the possibilities for which I^'^^ 
can have 0,1 or 2 elements. 

5.2. Special subspaces. For existence of non-trivial second order operators in 
Dy, it is necessary that / can be split into 2 chains. The form of T2 above is 
under the assumption that / can be split into two chains in exactly one way. 
There are cases, where / can be split in more than one way. As an example 
consider / = {0, 1, 2, ... , 98, 100}. Then / = {0, 1, 2, ... , 98} U {100} or / = 
{0, 2, 4, . . . , 98, 100} U {1, 3, 5, . . . , 97}. Consequently the space V admits a more 
general second-order operator. Here we describe all special cases, which fall into 4 
groups: 

Case A. The dimension oiV is 3, so I ^ {0,m,m + I}, I m. Then 

T2 = aix^+"'{xd ~m){xd-l-m) + a2X^+^d{xd - I - m) 
+ azx^ixd — m){xd — I ~ m) + a^x^d^ + a^xd + 

T2 gets a certain extra terms in the following cases: 

(a) m = 1, / > 2. Extra term aYd{xd — I — m). 

(b) m = 1, / = 2. Extra terms ard{xd — 3) + asd"^. 

(c) 1 = 1 (hence m > 1). Extra term o-id{xd — m). 

Case B. The dimension of V is 4, and / is "symmetric" , i.e. / = {0, m, m+l, 2m+l}, 
I ^ m. Then 

T2 = aix'+'^ixd - 2m - l){xd - I - m) + a2x"'ixd - 2m - l){xd - m) 
-|- a^x^d^ + a^xd + 0:5 

T2 gets an extra term only if 

(a) m = 1. Extra term aed{xd — m — l). 
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Case C. / has one runner ahead at distance 2, i.e. / = {0, 1, 2, . . . , n — 2, n}. Then 
T2 takes the form 

T2 = aix^{xd — n){xd — n + 3) + a'2x{xd — n){xd — n + 2) 
+ a^x^d^ + Uixd + 
+ aQd{xd — n) + a^d^ 

Case D. / has one runner left behind at distance 2, i.e. / = {0, 2, 3, 4, . . . , n}. Then 
T2 is of the form 

T2 = aix^{xd — n){xd — n + 1) + a2X^d{xd — n) 
+ Oizx^d"^ + Oi^xd + a5 
+ af,d{xd - 2) 

These special cases exhaust the list of all exceptional cases, which do not belong 
to the general classification given above. Let us prove this. It is clear that we 
can assume that the dimension is 5 or more, since the dimensions 3 and 4 are 
considered above. We know that / can be split into two chains of step, say, m, so 
*3 = «i — m, is = «3 — rn, . . . , i2r-\ and = i2—m, ie = 14— m, . . . , i2s- Suppose 
that / can also be split in two chains of step I. We may assume that I > m. Now 
we distinguish three cases: 

(1) r = 1, so 43 is not present. 

We have either /^'^ D {11,12} or j''^ D Therefore we consider two 

subcases: 

a. D {«i,«2}, so I ^ ii — 12- Since | I*^') |= 2, we have ^ so 
ii — in = I. But then — ig ^ and therefore 12 — ie = 2m = I. 
This is not allowed, since i\—ii = 2m implies that ii — 12 = fn. 

b. J*^') D {11,14}, so zi — ^2 = I- But also ^ and this is only 
possible \l 12 — i& = I = 2m. Therefore this configuration leads to 
case C. 

(2) s = 1, so 14 is not present. 

Clearly D {i\,iz}. Moreover i2 < 13 since otherwise also i2 G But 
then is ^ J^'^ implies that 15 + 1 = ii, and hence I = 2m. Again using 
i2 ^ ^^'^ implies that 12 + 1 = i2r-i- Therefore this configuration leads to 
case D. 

(3) r > 1 and s > 1, so is and 14 are both present. 

Always /^'^ D {ii, is}, and necessarily i2 < is, since otherwise also i2 € 
Again we have two subcases: 

a. is is present. We need is + Z = ii, sol = 2m. Like in case 2 above, 
it follows that 12 = i2r-i — 2m,. But then (4 G 

b. ie is present, but not is. We need either ii — 12 = I ov — 12 = I. 

If ii — i2 = I, then ig + / = i2 is the only possibility, so I = 2m, and 
hence is — i2 = m. This is forbidden because then we have one chain. 
If is — 12 = I, then i4 + I ^ {ii,is,i2}, so i4 G I'^^K So this is also 
impossible. 
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6. Conclusion 

In the previous sections we considered the case that the powers of x arc natural 
numbers. From algebraic point of view, this is not a crucial restriction. One could 
take the powers to be all integer, rational, real or complex numbers, or any other 
abelian subgroup of C. (If one takes C, one has to define a suitable total ordering on 
C to be able to consider positive and negative). This has as a main advantage that 
such algebras of "polynomials" admit the isomorphism x ^ \/x. Moreover, the 
map x'^ I— > a;'^+' is a linear isomorphism; it is the gauge transformation discussed 
before. 

An even more special case is that the set of allowed powers form a field. In this 
case, the change of basis x ^ .t™ is an isomorphism. 

We shortly discuss the case that the powers of x are real numbers. The algebra 
of "generalized polynomials", which we denote (suggestively) by C[x°'] is as a linear 
space 

and the multiplication is given by 

xl^ ■ x'^ = x'^+'^ /3,7GM 

Let D denote the algebra of differential operators with coefBcients in C[a;"] now. 
The degree of T e 23 can be any real number. For T e 2), deg(T) = m we have 

k 

T = ^c,x'+"'d\ 

i=0 

with Cj G C, and Cfc ^ 0. Lemma 3.1 for this operator T now looks like: 

Lemma 6.1. There exist numbers ai, . . . , G C such that 

T = CkX^{xd — ai){xd — 02) ■ ■ ■ {xd — a^). 

In particular, the case that m < dissappears, or, better, factorizes: 

d' = x-'{xd - (s - l)){xd - (s - 2)) • • • {xd - l)xd {s > 0) 

We denote I = {ii,i2, ■ ■ ■ , in}, where ij G M, j = 1,. . . , n, and define for m G M 

jt™) = {i^i\i + rn^I}. 

Let V = . . . and let Dy denote the algebra of operators in S) that 

leave V invariant. One can easily check that Theorem 3.2 holds literally. Due to 
the remarks at the beginning of this section, we have the following. 

Proposition 6.2. 

(1) The gauge transformation T 1— > x'Tx"' gives an isomorphism between 
andDw, with V = {x'\x'\... ,x'^) andW = {x'^-\x'^-\ . . . ,x'"-'). 

(2) The change of basis x' = x"^ induces an isomorphism between Dy and T)w, 
where V = {x^\x^\... ,x'^) andW = (a;'i/™,a;'^/'", . . . ,x^^/'^). 
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Proof. 

(1) Part (1) is obvious. 

(2) By the change of basis x' = x"^ we have that 

xd = mx'd', d' = 

ax' 

Hence the operator T = T-{xd—ai){xd—a2) ■ • ■ (xd—ak), where deg(T') = I 
is mapped to T' = T' ■ {mx'd' — ai){mx'd' — 0:2) • • • [mx'd' — au) with 
deg(T') = l/m. This leads directly to the proof of statement (2) of the 
Proposition. □ 

The classification of second order operators T2 G S)y is similar as before. A dif- 
ference is that l/^™)] = j/^""*)!, so that there is always a sort of symmetry in T2. 
Here we give one example, the generic 2-chain case. Suppose |/(™)| = 2, so 7 has 
the structure 7 = {ii,Z3, . . . ,i2r-i} U {i2,*4, • • • ,i2s\ with 

h = ii - m, = is - m, ... , i2r-i = hr-a - m 

and 

^4 = ^2 - m, if, = ii - m, . . . , i2s = «2.s-2 - rn 
If for alll > 0,1 ^ m there holds |/'-'^ | > 2, then the most general 2"'' order operator 
T2 S Dv takes the form: 

T2 = aix"^{xd — ii){xd — 12) + a2X^ + a^xd + a4, + a5X~"^{xd — i2r-i){xd — i2s) 

Note the symmetry in T2: there are as many terms of degree m as terms of degree 
— m. 

Final remark. In the case that V = {1 ), the algebra Dy is 

essentially generated by {d, 2xd — n, x'^d — nx}, see [4]. These 3 elements form the 
Lie algebra s^2(C), and hence Sy is essentially some representation of the universal 
enveloping algebra of 5/2 (C). We were not able to find a similar structure in J)y 
for general V, even for the simplest case V = {1, x, .t'^). Particularly, one can show 
that for the space V = {l,x,x^) the algebra Dy is an infinite-dimensional, finite- 
generated algebra. It is defined by 11 generators, which are differential operators of 
first-, second- and third order, and the commutator of any two of them is expressible 
as an ordered cubic polynomial in these generators. 
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